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We calculate the frequencies and damping rates of the 
low-lying collective modes of a Bose-Einstein condensed 
gas at nonzero temperature. We use a complex nonlinear 
Schrodinger equation to determine the dynamics of the con- 
densate atoms, and couple it to a Boltzmann equation for 
the noncondensate atoms. In this manner we take into ac- 
count both collisions between noncondensate-noncondensate 
and condensate-noncondensate atoms. We solve the linear 
response of these equations, using a time-dependent gaus- 
sian trial function for the condensate wave function and a 
truncated power expansion for the deviation function of the 
thermal cloud. As a result, our calculation turns out to 
be characterized by two dimensionless parameters propor- 
tional to the noncondensate-noncondensate and condensate- 
noncondensate mean collision times. We find in general quite 
good agreement with experiment, both for the frequencies and 
damping of the collective modes. 

PACS numbers: 03.75.Fi, 03.65.Db, 05.30.Jp, 32.80.Pj 



I. INTRODUCTION 

The nonzero-temperature dependence of the frequen- 
cies and damping of collective modes in trapped atomic 
Bose gases has been investigated extensively both ex- 
perimentally and theoretically [p|-p4||. Above the 
Bose-Einstein transition temperature the lowest-lying 
collective modes have been calculated in the collisionless 
regime using a collisionless Boltzmann or Landau- Vlasov 
equation Jl7| ] , and in the hydrodynamic regime either by 
using the conservation laws of hydrodynamics Q or by 
using the appropriate quantum kinetic equation fl. In 
the collisionless regime the frequency of the collective 
mode is large compared to the mean collision frequency, 
in contrast to the situation in the hydrodynamic regime. 
Several papers have also studied the regime intermedi- 
ate between the collisionless and hydrodynamic regimes 
by taking into account interatomic collisions fl7| , |23| , p4| , 
which turn out to mostly lead to mode damping. These 
papers have indicated that the experiments were per- 
formed under conditions intermediate between collision- 
less and hydrodynamic. (This is particularly true for the 
experiments performed by Mewes et al. Q and Stamper- 
Kurn et al. 0]). For temperatures far below the transi- 
tion temperature the collisionless modes can be described 
accurately by the time-dependent Gross-Pitaevskii equa- 
tion pa]. At higher temperatures the noncondensate 



fraction becomes substantial and the modes of the con- 
densate are now coupled to those of the thermal cloud. 
In Refs. [pp,|ri"| the temperature dependence of the mode 
frequencies has been calculated by employing the Popov 
approximation to include the static mean-field effects of 
the noncondensate atoms in the Gross-Pitaevskii equa- 
tion. This was improved later by taking into account 
also the dynamics of the thermal cloud using a RPA ap- 
proximation ] [Tof or the collisionless Boltzmann equation 
pl[ . For the hydrodynamic regime a two-fluid model has 
been developed in Refs. fl^ , [l5i. M oreover, the theory of 
Zaremba, Griffin, and Nikuni ]14| was improved later by 
the same authors to include collisions between conden- 
sate and noncondensate atoms |lfi| |. 

Using a combination of two previous papers, namely 
Refs. 21 1 and p4| , we aim in this paper at interpolat- 
ing between the collisionless and hydrodynamic regimes 
for experiments below the transition temperature. The 
authors of Ref . |p3 have already performed such an inter- 
polation for a gas above the critical temperature by using 
a Boltzmann equation with a relaxation time approxima- 
tion. In contrast, the authors of Ref. use a collision- 
less Boltzmann equation for the thermal cloud coupled to 
a time-dependent nonlinear Schrodinger equation for the 
condensate to consider the collisionless dynamics below 
the critical temperature. Here, we thus combine these 
two approaches by using again the Boltzmann equation 
and time-dependent nonlinear Schrodinger equation, but 
now including also the effects of interatomic collisions in 
the manner as put forward in Refs. @ and ||. This 
implies that we have to add two collision terms to the 
Boltzmann equation. The first collision term represents 
collisions between two noncondensate atoms and the sec- 
ond describes collisions between a condensate and a non- 
condensate atom. We use for both these collision terms 
a relaxation time approximation, since this approxima- 
tion leads above the transition temperature to a good 
agreement with microscopic calculations as well as with 
experimental data p3]. Furthermore, for consistency rea- 
sons we also have to include a damping term in the time- 
dependent non-linear Schrodinger equation, which is due 
to collisions between the condensate and noncondensate 
atoms. As a result our calculation will essentially be 
characterized by two parameters, namely T22 and T12, de- 
noting the mean collision time for collisions between the 
noncondensate atoms and between condensate and non- 
condensate atoms, respectively. Note that we use here 
the same notation as in Ref. |1^] . This will allow us to in- 
vestigate the collisionless and hydrodynamic limits with 
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respect to both t 2 2 and T12, and enable us to fit our re- 
sults for frequencies and damping with the experimental 
data. 

To solve the complicated nonlinear dynamics of the 
gas, we employ a gaussian trial function for the wave 
function of the condensate atoms with three complex 
time-dependent variational parameters. For the thermal 
cloud we use a distribution function that incorporates 
deviations from the Bose-Einstein distribution function. 
The deviation function is a truncated power expansion 
in the momenta and coordinates of the atoms. Our sys- 
tem of coupled equations is then assembled from the 
Euler-Lagrange equations for the variational parameters 
in the condensate wavefunction and from all the linear 
and quadratic moments of the Boltzmann equation for 
the distribution function of the noncondensate atoms. 
The solution of the linearized version of this system of 
coupled equations results in a dispersion relation that 
gives the frequencies and damping of the coupled modes 
in terms of T22 and t\%. Although we restrict our cal- 
culation to the axially symmetric traps used in the ex- 
periments, the generalization to fully anisotropic traps is 
straightforward. 

The rest of this paper is organized as follows: In the 
next section we present the theoretical details of our cal- 
culation. Specifically, we write down the time-dependent 
nonlinear Schrodinger equation for the condensate wave 
function and the linearized Boltzmann equation for the 
thermal cloud, and show how to treat the two collision 



terms in a relaxation time approximation. In section II A 
we present our trial functions for the condensate wave 
function and noncondensate distribution function and in 



section [I B we discuss how we obtain the equilibrium 



state of the gas around which we have to expand to find 
the collective modes. In section III we present the result- 



ing dispersion relation and discuss its collisionlcss and 
hydrodynamic limits. We then compare our results with 
experiment. We end in section ^ by summing up our 
main conclusions. 



II. COUPLED DYNAMICS OF THE 
CONDENSATE AND THE THERMAL CLOUD 

It has been shown that the coupled time-dependent 
nonlinear Schrodinger equation and the quantum Boltz- 
mann equation can both be derived starting from the 
equation of motion [R[lTJlHj23] 



lh —&r 



h 2 v 2 



+ V ext (r) 



2m 



(1) 



for the Heisenberg field operator ip(r,t). Here, g is the 
effective two-body interaction, which is given in terms 
of the scattering length a and the atomic mass m as 



g = 4irah 2 /m. The external potential we take here is 
a harmonic potential of the general form, 



V ext (r) = -ro(^ + w 2 V + ^3 A 



(2) 



where u>i is the characteristic frequency of the trap in the 
ith direction. To obtain the time-dependent nonlinear 
Schrodinger equation we write i/i(r,t) as 



^(r ) i)=*(r,t)+^ / (r,t), 



(3) 



where $(r, t) is the appropriate nonequilibrium expecta- 
tion value of V?(r, t) and the operator ip'(r, t) describes the 
noncondensate atoms. The number density of condensate 
atoms n c (r,t) is related to $(r, t) by n c (r, t) = |<&(r, t)\ 2 
while the number density of the noncondensate atoms n' 

equals (r, i)ifi'(r, t)Y where the brackets represent 
the averaging over the nonequilibrium density matrix. 
Substituting this form for ip( r ->t) m the equation of mo- 
tion and averaging, we obtain 



d $(r,t) 
dt 



h 2 V< 

2m 



+ V ext (r) + gn c (v,t) + 2gn'(r,t) 



-iR(r,t) J$(v,t), (4) 
where R(r, t) is given by 

fl(r,t) = - ff ((^(r,f$'(r,t))**(r,t) 

+ ($>\t, t)V>(r, t$'(r, i))) /$(r, t). (5) 

Next, the quantum Boltzmann equation can be derived 
by writing an equation of motion for the distribution 
function of the noncondensate atoms /(p,r,i). This is 
usually done by writing /(p, r, t) as a Wigner transform 

of (r', t)ip'(r, t)^ and then determining the time evo- 
lution of /(p, r, t) from the equation of motion in Eq. (|l|) 
p8| , p2 J26 1 . In the Hartree-Fock approximation the result- 
ing Boltzmann equation takes the form 
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f = C 22 [f] + C 12 [f], 



(6) 



where C22 is the contribution to the rate of change of 
/ due to collisions between noncondensate atoms, while 
C12 is the contribution due to collisions between the 
condensate and the noncondensate atoms. The energy 
E(p, r, t) of the noncondensate atoms is in this approxi- 
mation given by |2l| l 

E(p, r, t) - p 2 /2m + V ext (r) + 2gn(r, t), (7) 

where n — n c + n' is the total density. 

In Ref. explicit forms of the two collision terms 
have been written down as follows 
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C 22 [f] = (2^)5^ / ^ j j dPi5 ( P + P2 - P3 - P4) 

x S(E + E 2 -E 3 - E A ) 

x [(1 + /)(1 + f 2 )hh - // 2 (1 + / 3 )(1 + U)\ , (8) 



and 



Cl2[f] 



dpi / dp 2 / dp 3 



2g 2 n c 
~ (2tt) 2 /i 4 
x 6(mv c + pi - p 2 - Ps) 
x 6(E C + E X -E 2 - E 3 ) 
x [8(p - pi) - S(p - p 2 ) - <5(p - p 3 )] 
x [(l + /i)/ 2 /3-/i(l + / 2 )(l + / 3 )]. 



(9) 



Here E c (r,t) and v c (r, t) are the local energy and ve- 
locity of the condensate atoms, Ei is a shorthand for 
E(pi,r,t), and similarly fi is a shorthand for /(pi,r, t). 
As is shown in Ref. |l8[ ], the conservation of the total 
number of atoms constraints the collision term C\ 2 in the 
Boltzmann equation to be related to the damping term 
R(r, t) in the complex nonlinear Schrodinger equation as 



2n c (r,t)R(r,t) 



dp 



■Ci 2 [/]. 



(10) 



The coupled equations given in Eqs. (|J), (||), and jlO| ) 
in principle fully determine the dynamics of the Bosc- 
Einstein condensed gas in the Hartree-Fock approxima- 
tion. Since the dynamics of the thermal cloud is ex- 
perimentally only important for temperatures which are 
larger or comparable to the mean-field interactions we 
believe that this is sufficiently accurate and we do not 
need to use the Popov approximation at first instance. 
Nevertheless, these equations are too complicated to be 
solved exactly and some approximation is called for. 



A. Trial Functions 

For the noncondensate atoms we start by linearizing 
the Boltzmann equation in small deviations of the distri- 
bution function around its equilibrium value, namely 



/(p, r, t) = fW (p, r) + /W (p, r) (1 + fW (p, r))^(p, r, t), 

(11) 

where f^°'(p,T) is the equilibrium distribution function 
given by 



/W (p, r) = {exp [(E(p, r) - fi/k B T\ - l}" 1 



(12) 



Here /z is the chemical potential and ks is Boltzmann's 
constant. As a result the linearized Boltzmann equation 
takes the form 



^ + V p £ • V r - V r E ■ V p 



tP = C 22 [iP} + C 12 Mj] 



It is shown in Ref. p4| that for an uncondensed gas a 
trial function for ip(p, r, t) of the form 



ip = Aix 2 + Bixp x + C'ipl + A 2 y 2 + B 2 yp y 
+ C 2 pl + A 3 z 2 + B 3 z Pz + C 3 p 2 z , 



(14) 



where Ai(t), Bi(t), and Cj(i) are nine time-dependent 
functions, is appropriate to describe the low- lying colli- 
sionless breathing modes with frequencies 2wj, some of 
which have been observed experimentally. We will there- 
fore assume the above expansion to be also reasonably 
accurate below the critical temperature. For the colli- 
sion integrals C22 and C\ 2 we use a relaxation time ap- 
proximation. In such an approximation one associates 
these collision integrals to mean relaxation times. How- 
ever, the approximate expressions for C 22 and C\ 2 should 
still take into account the conservation laws associated 
with the collision processes exactly. In the case of C 22 , 
which represents collisions between the noncondensate 
atoms, the number of atoms in the thermal cloud, their 
total momentum and their total kinetic energy should 
all be conserved. In our ansatz for ifj, terms like x 2 and 
xp x should therefore not be affected by such collisions 
since they correspond to two collision invariants, namely 
the number of atoms and the total momentum in the x 
direction, respectively. On the other hand, terms like 
p 2 , will be affected by collisions since it is only the sum 
p 2 = p 2 . +Py+ p 2 z which is conserved during the collision 
process. Therefore, we can write the following expression 
for the linear operator C 22 [|2C)|,p4[| 



C 22 [ipi] = 



1 

T-22 







otherwise 



(15) 



where t 22 is a mean collision time for the noncondensate- 
noncondensate atomic collisions. Note that in this ex- 
pression C 22 [p 2 ) = 0, ensuring the conservation of the 
total kinetic energy. 

For C12, which represents collisions between the con- 
densate and the noncondensate atoms, the number of 
atoms is not conserved, since the collision process in- 
volves transport of atoms back and forth from the con- 
densate into the thermal cloud. This statement means 
mathematically that the zeroth moment of C\ 2 , i.e., 
J dpC\ 2 , does not vanish in contrast to the case of C22 
where J dpC 22 — 0. As a first attempt to associate 
with Ci 2 a mean collision time T\ 2 for noncondensate- 
condensate atomic collisions, we may write C\ 2 cx \ jr\ 2 . 
However, we observe from Eq. that C\ 2 cx n c (r). 
This dependence on n c (r) requires us to assign to r\ 2 a 
position dependence that follows from ro c (r) as 



1 _ 1 n c (r) 
ri 2 (r) n 2 (0)n c (0)' 



(16) 



(13) 



where ti 2 (0) and n c (0) are the noncondensate- 
condensate mean atomic collision time and condensate 



3 



density at the center of the trap, respectively. Generally, 
Ti2(0) may, just like T22, also have temperature depen- 
dence. By multiplying and dividing l/ri2(r) by kg we 
can thus take C\ 2 as obeying 



C12 oc a 



2gn c {v) 



where we introduced the dimensionless constant 

h 



2 3 n c (0)r(0) 



(17) 



(18) 



This form for C12 is convenient for a reason that be- 
comes clear lateron in this section. It should be empha- 
sized here that the dimensionless parameter a is the sec- 
ond free parameter in our phenomenological calculation. 
The first dimensionless parameter being 1/Qt 22 where 
Co = ^loioj 2 lo 3 . Therefore, the collisionlcss and hydrody- 
namic regimes, as well as the intermediate regime, will 
be described by only these two parameters. Finally, we 
make the relation between C12 and a precise by taking 



2gn c (r) 




4>i 



2 2 2 
Px,Py,Pz Z 

tpi = x 2 ,y 2 ,z 2 
otherwise 



(19) 



This completes our description of the treatment of the 
thermal cloud. Next we have to consider the condensate. 

It is known that a time-dependent gaussian ansatz for 
the condensate wave function gives the correct frequen- 
cies of the lowest modes at zero temperature p7j-gj| ■ Fur- 
thermore, the gaussian ansatz has also been used in Ref. 
pl[ for the whole temperature range below the transition 
temperature and leads to rather good agreement with ex- 
periment. Therefore, we employ here again a gaussian 
ansatz for the wave function of the condensate. It has 
the following form 



$(r,t) 



8N c bi r b 2r b 3r 



exp [-(hx 2 + b 2 y 2 + b 3 z 2 )], 

(20) 



where 61, b 2 , and b 3 are complex time-dependent vari- 
ational parameters and b± r , b 2r , and b 3r are their real 
parts, respectively. Similarly, we denote the imaginary 
parts, which will appear lateron, as bu, b 2i , and 63^. 
The prefactor of $(r, t) guarantees its normalization 
J dr|<&(r, t)\ 2 to be equal to the number of condensate 
atoms N c . 

To obtain the first set of our coupled equations of mo- 
tion, we start by writing down the energy functional that 
corresponds to the nonlinear Schrodinger equation in Eq. 
(g), namely 



E c [$] = Jdr 



h 2 |V$| 
2m 



1 



+ V ext n c + -g(n c 



4n c n') 



— iRn r 



(21) 



with n c = \$\ 2 . We evaluate this energy functional using 
the expression for <& from Eq. (|20|), and by noting that 
the density of the thermal cloud 



dp 

W) 1 



/(P,r,t), 



(22) 



and the damping term i?(r, t) can be calculated using 
Eqs. jlO| ) and (|l9|). We notice here that —iR(r,t)n c 
turns out to be proportional to the Hartree-Fock inter- 
action term in Eq. (|2l]), with a proportionality constant 
equals to —ia. This explains the reason for writing C\ 2 
as in Eq. (fl7|). The equations of motion for the con- 
densate dynamics are now the Euler-Lagrange equations 
resulting from varying the lagrangian 



-ih 



>h- j 

at 



<4*. 

dt 



em 



(23) 



with respect to the 6 variational parameters bk r and bki ■ 
The second set involves the equations of motion for 
the constants Ai, Bi, and Cf. It is obtained from taking 
appropriate moments of the Boltzmann equation in Eq. 
(|l3|). In detail, the moments are calculated by multiply- 
ing this equation by /W(p, r) and the various compo- 
nents of ip in Eq. (|l4|), and then integrating over p and 
r. This results in 9 equations of motion. In combina- 
tion with the previous 6 ones, we thus have 15 coupled 
equations of motion. The coupling is provided on the 
one hand by the mean-field interaction 2gn c n' and the 
imaginary damping term iRn c in the condensate energy 
functional E c in Eq. (|2l]), and on the other hand by the 
contribution 2gn c to the Hartree-Fock energy E in Eq. 
(0) and the C\ 2 collision term in the Boltzmann equation. 



B. Equilibrium 

Finally, it is important to note that linearization of 
the equations of motion of the condensate is obligatory 
to be consistent with the equations of motion of the non- 
condensate part which are already linearized. Therefore, 
we need to calculate the equilibrium state of the gas. 
In principle, to obtain equilibrium properties we should 
minimize the free energy F — E to t — TS, where E to t is 
the total energy and S is the entropy of the gas, with 
respect to some variational parameters that character- 
ize the widths of the condensate and the thermal cloud. 
However, a simplified estimate of the total energy and 
the entropy contribution to the free energy above the 
transition temperature shows that for the experimental 
conditions of interest the former is dominant over the 
later. This simplified estimate can be made by using a 
one-parameter gaussian ansatz for the thermal cloud den- 
sity, namely n' oc exp(— x 2 /R 2 ), where R is the radius of 
the cloud. It turns out that in the Thomas-Fermi limit 
the entropy contribution to the free energy is of order 
/csTln7 3 / 5 , where 7 = N'a/a, N' is the total number 
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of atoms in the thermal cloud, and a = y'fi/mw is the 
harmonic oscillator length, while the total energy con- 
tribution is of order ksTj 2 / 5 . Therefore, if 7 3> 1 the 
entropy contribution is much less than the total energy 
contribution. 

As we shall see in the next section, the experi- 
ments were performed with a temperature-dependent to- 
tal number of atoms ranging from about 6000 atoms at 
zero temperature to approximately 40000 atoms at the 
transition temperature. If we extend this simple esti- 
mate to temperatures below the transition temperature 
it turns out that the condition 7 3> 1 is satisfied only 
for high temperatures but not for low temperatures since 
the number of thermal atoms becomes small. However, 
this condition is then no longer important since the free 
energy of the thermal cloud at such temperatures is small 
compared to the condensate energy. Therefore, we will 
in the following only minimize the total energy E to t with 
respect to the variational parameters of the condensate 
and the thermal cloud. To be able to calculate E tot we as- 
sume that the distribution function of the noncondensate 
atoms has the same form of that of a noninteracting gas 
but with varying spatial widths that effectively take into 
account the mean-field effects of both the noncondensate 
and condensate atoms. This effective distribution func- 
tion is written as 

^ff-[ eXP { k B T + i? 1 2 + i? 2 2 + i? 3 2 J 

(24) 

where R\, R2, and R3 are the widths of the nonconden- 
sate cloud in the three directions. The total energy is 
now a function of 6 variational parameters: 

Etot — Etot {b^iJ , &2r > ^3r > j -^2 > -R3 ) j (25) 

where b^J , b^J , and bfj are the equilibrium values of b\ r , 
b 2r , and 6 3r , respectively. The equilibrium is obtained by 
minimizing this energy with respect to these variational 
parameters. 

The results of such a minimization will be shown in 
the next section, where we present the dispersion relation 
that results from solving the above-described system of 
15 linearized coupled equations of motion. We discuss 
also the coUisionless and hydrodynamic limits of these 
results. 



m the projection of the angular momentum of the mode 
along the axis of symmetry of the trap . We focus in 
this paper on the two lowest-lying m = and m = 2 
modes observed experimentally. It turns out that for 
the experimentally relevant temperature range, the in- 
phase m = 0, and m = 2 modes correspond mostly 
to oscillations of the thermal cloud, whereas the out-of- 
phase modes are mostly condensate oscillations. There- 
fore, we shall often refer in this paper to the in-phase 
modes as the thermal cloud or noncondensate modes 
and the out-of-phase modes as the condensate modes. 
Although our calculation provides results for another, 
higher-lying, monopole mode, we shall not discuss it 
further here. Moreover, we emphasize that throughout 
the following we perform our calculations for parameters 
taken from the experiments of Jin et al. Q , i.e., with 87 Rb 
atoms in an axially-symmetric trap with anisotropy ratio 
uj 3 /lji — w 3 /oj 2 = 

In that particular experiment the measurements were 
performed with a temperature-dependent total number 
of atoms as a result of the loss of atoms during evapo- 
rative cooling. The total number of atoms N tot , as well 
as the number of condensate atoms -/V c , are measured in 
the temperature range T/T BEC ^ 0.48 to T/T BEC ^ 1.0, 
where Tbec is the Bose-Einstein transition temperature. 
These measurements can be easily fitted with polyno- 
mials in T/Tbec as shown by the dashed line in Fig. 
|l|. An extrapolation of such a fit to temperatures below 
0A8Tbec leads, however, to nonphysical situations. In 
particular, the two curves for N tot and N c cross at least 
once before reaching zero temperature. We overcome this 
problem by fitting only N to t with the experimental data, 
and then using N c —N to t(^ — (T/Tbec) 2 ) m analogy with 
the ideal gas relation N c =N tot (l - (T/T B ec) 3 ) %■ Fig. 
|l| shows with the solid line the results of this slightly 
less accurate fit. It should be noted that using the ideal 
gas relation will grossly overestimate all the experimen- 
tal points of N c for T/Tbec > 0.7. It turns out that the 
difference between the calculated frequencies using the 
less accurate fit from those calculated using the best fit 
are much smaller than the uncertainties in the measured 
frequencies for the experimental range of temperatures. 
This is shown in Fig. ^[ Since we want to show also 
results for the complete temperature interval from zero 
to T B eCi we employ from now on always the former fit, 
i.e., the solid lines in Fig. pi. 



III. THE DISPERSION RELATION: 
FREQUENCIES AND DAMPING RATES 

Our calculation accounts, in fully anisotropic traps, 
for 9 modes of the gas. In axially symmetric traps this 
number reduces to 6 modes. These are the in-phase and 
out-of-phase combinations of the two monopole (m = 0) 
modes and one quadrupole (m = 2) mode of both the 
condensate and the thermal cloud. Here we denote with 
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FIG. 1. The measured total number of atoms Ntot 
(squares) and the number of condensate atoms N c (circles). 
The dashed lines are a polynomial fit in powers of T/Tbec to 
the experimental data . The solid lines represent a fit to Not 
only and using for N c the relation iV c = N tot (1 - (T/T BEC ) 2 ). 
The solid line fit of Ntot is shifted slightly downwards in order 
to get better agreement with the data for N c . The measured 
points were taken from Ref. [2]. 




BEC 



FIG. 2. Collisionless frequencies for a = 0, i.e., with no 
collisions between condensate and noncondensate atoms. The 
solid lines correspond to the approximate fit to the total and 
condensate number of atoms as shown by the solid lines in 
Fig. The dashed lines correspond to our best fit in the 
temperature interval as is given by the dashed lines in Fig. [IJ 
The points are the experimental measurements for the m — 
and m = 2 condensate modes [2] . The labels of the theoretical 
curves in terms of the values of m are shown in Fig. ^ 

Before starting with calculating the frequencies and 
damping rates of the collective modes, we show in Fig. 



H the result of the minimization of the total energy that 
is required to obtain the equilibrium conditions of the 
gas, as described in the previous section. We plot the 
equilibrium widths of both condensate and nonconden- 
sate clouds as a function of temperature. We notice from 
this figure that at zero temperature and at the transi- 
tion temperature we obtain rather good estimates for 
the widths of the condensate and thermal clouds. At 
zero temperature, where the thermal cloud is absent, we 
observe that the radial width of the condensate equals 
approximately 2ai, a result that is consistent with the 
familiar property that due to the mean-field interactions 
the equilibrium condensate width is larger than the ideal 
gas result a\. The axial width, which is roughly 1/A 1 / 2 of 
the radial width, is suppressed due to the anisotropy of 
the trap. We can also see in this figure that at the tran- 
sition temperature the radial width of the condensate is 
slightly larger than the ideal gas one. This slight expan- 
sion is caused by the presence of the thermal cloud, whose 
mean-field interaction has the effect of slightly reducing 
the spring constants of the effective trapping potential. 
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FIG. 3. The equilibrium widths of the condensate 
and noncondensate clouds in the axially symmetric trap 
of Ref. [2]. In reference to the text, the radial width 
of the noncondensate cloud is Ri(= R2), and the axial 
width is 7?3. For the condensate cloud, the radial width is 

1/2/ (Oh^h- 1 / 2 ^ nnr] the avial wiH+.Vi is ( ')h < - '> \~ 1 / 2 



(2b ( °J)-' rz (= (2b%>)- L/z ) and the axial width is (26^' r 
as can be seen from Eq. (p5|). All lengths here are scaled to 
the harmonic oscillator length a\ = (h/mcui) 1 ^ 2 . 



Returning to the problem of the collective mode fre- 
quencies and damping, we start by neglecting collisions 
between the condensate and noncondensate atoms, which 
means that a = 0. In this case the dispersion relation 
turns out to have the following general structure 



iPf>UJTZ2 



P, 



(0: 



H 



0. 



(26) 



where Pq\ P\ - , and P^' are 6th order polynomials in 



,(0) 
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u! 2 which can all be factorized as 



p(0) _ 



IK' 

k=l 



P^ = 



k=l 



J Ck) 



J Ik) 



and 



P. 



(0) 



H 



m 

k=l 



J Hk) 



(27) 



(28) 



(29) 



Here, u>ck, ^>Hk, and ujik, k = 1,2,3, are temperature- 
dependent collisionless, intermediate, and hydrodynamic 

frequencies, respectively. The superscripts in Pq\ Pj , 

and P^ indicate the value of a. 

The general structure of this dispersion relation agrees 
with the result of Refs. above the transition tem- 

perature, where this dispersion relation was studied in 
detail. Particularly, it was shown that damping rates 
calculated using this relation agree in order of magnitude 
with experiments and numerical calculations. 

The collisionless regime is defined by the condition 
WT22 3> 1. Therefore, we find from Eq. (|2^) that the 
collisionless frequencies are u>ck- In the hydrodynamic 
limit UIT22 <C 1, and the hydrodynamic frequencies are 
^Hk- Note that we use the word hydrodynamic here to 
denote that the thermal cloud is in the hydrodynamic 
regime though a = and there are therefore no colli- 
sions between condensate and noncondensate atoms. In 
a sense this regime is thus precisely the limit discussed 
by Nikuni, Zaremba, and Griffin Using the experi- 
mental parameters we calculate these frequencies for the 
whole temperature range below Tbec- In Fig- § we 
present the results of this calculation for the m — and 
m = 2 modes together with the experimental data. In 
this figure the collisionless curves agree with those of Bi- 
jlsma and Stoof Q for T /Tbec > 0.2, which is not 
surprising since we use similar ansatz functions. The dis- 
crepancy for T/Tbec < 0.2 is due to the different ways 
of treating the equilibrium state. 

Above the transition temperature, the explanation of 
the mode damping was based on the fact that the mea- 
sured frequencies were less than the theoretical collsion- 
less frequencies. This shift in frequency is then inter- 
preted to be due to the fact that the system is shifted 
from the collisionless regime towards the hydrodynamic 
regime. Collisional damping associated with this shift 
was calculated and compared to the experimental damp- 
ing, ft is clear from Fig. [| that below the transition 
temperature this kind of explanation is not possible, since 
most of the experimental points are not located between 
the collisionless and hydrodynamic curves. 



- 1.7 




-r 1-5 



FIG. 4. Collisionless (solid curves) and hydrodynamic 
(dashed curves) collective mode frequencies for a — 0. The 
curves are labeled in terms of the values of m and the phase 
difference between the condensate and thermal oscillations 9: 
(a) m = 0, = 0, (b) m = 2, 6 = 0, (c) m = 0, 6 = vr, (d) 
m = 2, 6 = 7T. 

To investigate the possibility that this damping might 
be due to the condensate-noncondensate atomic colli- 
sions, we now study the effect of these collisions by taking 
a to be nonzero. For nonzero a the dispersion relation 
takes the general form 



\k=0 



(wr 22 )' 



\k=0 



(fe) a fe 



UJT 2 2 



£4 fe v = o. 



(30) 



k=0 



We note that the k = term in this expression corre- 
sponds to Eq. (p6|), where a was indeed taken to be 
zero. Frequencies and damping rates can now be ob- 
tained by writing oj in terms of a real and imaginary 
parts, uj — uj r + iuii, and then inserting this expression in 
the last equation. In general, one obtains expressions for 
ui r and uii in terms of UJT22 and a. 

We start our discussion of the frequencies and damp- 
ing rates for T/Tbec — 0.79, since at this temperature 
the experiment provides data for both the condensate 
and noncondensate oscillations. We present the results 
of our calculation for the m = and m = 2 modes of 
the thermal cloud in Fig. ^, where we have eliminated 
UJT22 and plotted u>i as a function of ui r for a = and 
a = 0.25. On the same plot we show the two experi- 
mental points corresponding to the m = and m = 2 
modes of the noncondensate oscillation. In this figure, 
the right end of the curves represent the collisionless 
regime (cZ>t 2 2 ^> 1) and the left end is the hydrodynamic 
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regime (0JT22 <C 1). Therefore, the location where the 
dashed curves intersect with the horizontal axis can also 
be seen in Fig. | at T = 0.79T B ec- For the experi- 
mental point corresponding to the monopole mode, the 
error bars do not intersect with the dashed curve (a = 0) , 
but they do with the solid curve (a = 0.25) indicating 
that a nonzero value of a is needed to account for the 
experimental damping. This indicates that condensate- 
noncondensate atomic collisions are mainly responsible 
for the observed damping of this mode. This is less clear 
for the quadrupole mode although the experimental data 
is certainly not inconsistent with this conclusion. 

Next we show in Fig. ^ at the same temperature both 
the condensate as well as the noncondensate to — and 
to = 2 modes for values of a ranging from to 0.25. 
The effect of a on the noncondensate modes is a slight 
shift and rotation of the curves but their general struc- 
ture is roughly preserved. This can be seen more clearly 
in Fig. [|. The effect of a on the condensate modes is 
much larger, it gives rise to a large upward shift for the 
whole curve. The effect of 0JT22 on these curves becomes 
smaller for larger values of a, which can be seen by noting 
that the radius of the small semicircles becomes smaller 
for larger values of a. Again, we see that a nonzero value 
of a provides sufficient damping to account for the ex- 
perimental observations. However, there is a discrepancy 
between our calculation and the experiment in the value 
of the frequency of the m = mode, since our calcula- 
tion predicts lo/uji « 1.68 ± 0.4, whereas the experimen- 
tal data give a frequency ps 2.05a->i. This discrepancy 
was also present in the calculations of Refs. |l0|,^l| and 
is discussed in Ref . . The authors of the last reference 
suggested that a possible reason for this discrepancy may 
be that the observed mode is in fact an in-phase to = 
mode (the upper solid curve in Fig. rather than an out- 
of-phase one. This suggestion was based on calculating 
the oscillator strenghts of the in-phase and out-of-phase 
to = modes. It turned out that it is indeed possible ex- 
perimentally to excite both modes simultaneously in the 
temperature range T w 0.25T B ec to T m 0.5T BE c- Al- 
though the calculation in Ref. J21] was performed in the 
collisionless limit, where there is no damping, we expect 
that with damping this argument remains qualitatively 
correct. 




FIG. 5. Damping rate u>i versus frequency uo r for noncon- 
densate oscillations, at T = .79Tbec for a — and a = 0.25. 
The two points correspond to the experimentally observed 
frequency and damping rate of the thermal cloud oscillation. 




FIG. 6. Damping rate u)i versus frequency u r of conden- 
sate and noncondensate oscillations at T = .79Tbec, for a 
changing from to 0.25 in steps of 0.025. The two points cor- 
respond to the experimental measurements of the condensate 
oscillation. 



For the rest of the experimental points we can, in 
principle, explain the data using a two-parameter fit, 
namely with a and uiT22- In fact, the collisions between 
atoms from the condensate with those from the ther- 
mal cloud are the main cause of damping in the con- 
densate modes, and we can even explain the data us- 
ing a one-parameter fit, namely only a, by assuming 
that the system is in the collisionless regime with respect 
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to <I>T22, i-e., C0T22 3> 1. We perform this calculation 
by using a function a(T /Tbec) that gives the best fit 
for the experimental damping rates. In Fig. [?] we plot 
the function a{T /Tbec) and the resulting temperature- 
dependent frequencies and damping rates. This figure 
contains the main results of this paper. First of all we 
notice that now, with a nonzero a, we obtain better 
agreement with the experimental data for the quadrupole 
mode frequencies than before where a was taken to be 
zero. This can be seen by comparing this figure with 
Fig. ||. We note also that, depending on the value of a, 
the mode frequency may shift upwards or downwards. 
By comparing Fig. |] with Fig. ^ for the mode fre- 
quencies one can clearly see that for temperatures be- 
low approximately 0.6Tbec the theoretical curves are 
shifted upwards, whereas for the higher temperatures the 
curves are shifted downwards, in the end giving rise to 
the good agreement of the quadrupole mode with ex- 
periment. Secondly, we notice that the same function 
a(T /Tbec) also gives a good fit for the damping rates of 
both the monopole and quadrupole modes. Finally, there 
is still a discrepancy between our predictions and the ex- 
perimental findings for the frequencies of the monopole 
mode at higher temperatures, but this may be resolved 
by a reinterpretation of the observed mode, as explained 
previously. 




0.001 



FIG. 7. Fitting all the experimental data on the 
mode frequencies and damping rates with a single temper- 
ature-dependent function a. 

Finally, we want to mention an interesting feature of 
the in-phase and out-of-phase m = modes which can 
already be seen in Fig. || at T/Tbec ~ 0.05. It is the 
familiar anticrossing tendency of these modes obtained 
also by Bijlsma and Stoof 21 and Zaremba et al. |l8| . 
It is more dramatic here than in these two works due to 
the presence of the new parameter a. This behavior can 
be seen in Fig. ^ where we see the two modes come close 
to each other and eventually intersect at some 'critical' 
values of a twice. 
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FIG. 8. Damping rate u)i versus frequency u) r of conden- 
sate and noncondensate oscillations at T = .79Tbec with a 
changing from 0.0005 to 0.05. The four curves in each plot cor- 
respond to the in-phase and out-of-phase m = and m = 2 
modes which are also shown in Fig. [| The present figure 
corresponds to the region in Fig. M where the in-phase and 
out-of-phase m = modes come close to each other. One 
can see here that the anticrossing feature of these two m — 
modes is not only present at a = [21] but also for nonzero 
values of a. 



IV. CONCLUSION 

We have extended the method of calculating collisional 
damping, used previously above the Bose-Einstein con- 
densation transition temperature p3|j24| ], to below the 
transition temperature. Furthermore, we have included 
the effect of noncondensate-condensate collisions. By 
comparing with experiment, we conclude that it is pre- 
sumably this collision process which is mainly responsible 
for the observed damping. Our theory provides a general 
dispersion relation that gives complex mode frequencies 







at a certain temperature as a function of two dimension- 
less parameters, namely 1/lotzi and a, that characterize 
the noncondensate-noncondensate and noncondensate- 
condensate atomic collisions, respectively. Our results 
for the fully collisionless frequencies agree with those of 
Rcf. for most of the temperature range below the 
transition temperature. The m = and m = 2 hydro- 
dynamic frequencies are, to the best of our knowledge, 
calculated here below the transition temperature for the 
first time. 

At present we have not carried out a microscopic 
calculation of the mean collision times T22 and T12 oc 
a -1 , which arc treated as phcnomcnological parame- 
ters here with the possibility to investigate the inter- 
mediate regime with respect to either of the two colli- 
sion processes. However, the homogeneous calculations 
of Zaremba, Nikuni, and Griffin ]l2fl indicate that the 
function a(T /Tbec) that we have obtained from a fit 
to the experimental data, has the correct order of mag- 
nitude and qualitatively also the correct temperature- 
dependence. Nevertheless such a microscopic calculation 
for the inhomogeneous experimental conditions of inter- 
est would be very desirable and is left for future work. 

We obtain a rather good agreement with the exper- 
imental results for the mode frequencies and damping 
rates apart from the discrepancy for the m = out-of- 
phase mode for T > O.ITbec- While the experiment 
shows an upward shift in the frequency with respect to 
its zero temperature value, we predict a downward shift. 
Resolution of this discrepancy seems to require more ac- 
curate experimental data, which will also provide decisive 
comparison with the theory presented here. 
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